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Abstract 

p ^ . The problem of matrix factorization motivated by diffraction or elasticity is 

^ [ studied. A powerful tool for analyzing its solutions is introduced, namely analytical 

continuation formulae are derived. Necessary condition for commutative factoriza- 
tion is found for a class of "balanced" matrices. Together with Moiseyev's method 
and Hurd's idea, this gives a description of the class of commutatively solvable ma- 
trices. As a result, a simple analytical procedure is described, providing an answer, 
whether a given matrix is commutatively factorizable or not. 
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1 Introduction 

{Sj i A matrix factorization problem (i.e. a problem of finding the factors Q + and Q~ providing 

the decomposition (CQ) for a known matrix G(k)) is usually motivated by an elasticity or 
a wave diffraction problem. Typically its formulation does not contain a requirement of 
commutative factorization. However, the possibility to perform a commutative factoriza- 
tion is usually studied carefully, since all known factorization methods are based on ideas 
connected with commutativity. 

There are two main methods for commutative factorization. The first one is based on 
the idea by Heins [T] who proposed to split the logarithm of matrix G additively. Further 
progress in this direction is connected with the names of Chebotarev [2] and Khrapkov [3] . 
In the paper by Khrapkov an explicit form of factorization for a certain class of matrices 
2x2 was found. 

However, sometimes Khrapkov's formula leads to exponential growth of the factors at 
infinity, and this is not acceptable for physical applications. To suppress this growth in 
some cases, a special technique was proposed by Daniele [I] and Hurd & Liineburg [5], 
who thus enlarged the class of explicitly factorizable matrices. A factorization to some 
other matrices can be obtained if the analytical continuation of the matrix is studied (see 
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Rawlins [6] and Hurd [7J). Matrices with dimension more than 2x2 have been investigated 
by Lukyanov [5]. 

Another technique of commutative matrix factorization is based on diagonalization 
of G and studying the eigenvalues as a multi-valued function. This approach has been 
developed by Cercignani [9], Liineburg [ID] . Moiseyev [11], Meister & Penzel [12] and 
Antipov & Silvestrov [13]. It is based mainly on the research of Zverovich [14], who 
described a method for solving scalar boundary-value problem on Riemann surfaces. 

It was noticed by Hurd [15] and Daniele [16] that all known methods to factorize a 
matrix are connected with the fact that the matrix commutes with a polynomial matrix. 

A question of finding a rational factor enabling to perform a commutative factorization 
for matrices 2x2 was discussed by Williams [17] ■ Also this question was stied in detais 
by Ehrhardt and Speck [T8] . 

Beside the exact methods, an interesting attempt to solve a matrix factorization prob- 
lem approximately was proposed by Abrahams [19] . 

Each of the works having either Kharapkov's or Moiseyev's context was dedicated 
to a distinct class of matrices, i.e. the starting point of such a work was a phrase like 
"consider a matrix having the following form . . . " . Sometimes, however, it is not easy 
to say whether a given matrix can be reduced to one of the known classes by algebraic 
manipulations. That is why, a separate and interesting issue is a classification of matrices 
with respect to factorization, i.e. finding a criteria, for example, of the possibility of 
commutative factorization. Such a theory is known for matrices 2x2 |2J, however the 
necessary condition for commutative factorization has been found in an "Ansatz" form, 
i.e. a matrix should have a specific representation including some entire (polynomial) 
matrices and some arbitrary functions as coefficients. 

A necessary condition for commutative factorization for matrices of arbitrary dimen- 
sion was studied by Jones in [20], however the author restricted himself to the case of 
matrices having distinct eigenvectors everywhere, while a typical matrix emerging in 
diffraction theory has branch points, i.e. it has distinct eigenvectors almost everywhere. 

The current paper is inspired mainly by the works of Antipov et.al. and the bright 
work of Hurd. The idea is to take an algebraic matrix G(k) and study the properties of 
the factors Q + and Q~ on their Riemann surfaces following a priori from the decomposi- 
tion ([1]), i.e. without constructing the factors explicitly. We found that the decomposition 
(CQ) taken together with the regularity conditions imposed on the factors define a unique 
Riemann surface of Q + or Q~ . Moreover, the values of, for instance, Q + taken on dif- 
ferent sheets are connected by simple algebraic relations. Thus, in Section 2 we obtain 
analytical continuation formulae. The relation obtained by Hurd is a particular case of 
such formulae. 

The problem for Q + reminds the functional problem for Abelian integrals, but while 
the value of an Abelian integral is increased by a constant when the argument is carried 
along a certain contour, the value of Q + is multiplied by a non-constant bypass matrix. 

In Section 3 we study the question of the possibility of commutative matrix factor- 
ization. To describe the matrices, for which the necessary condition for commutative 
factorization is fulfilled, we introduce the class of branch- commutative matrices G(k), i.e. 
the matrices whose values on different sheets of their Riemann surface commute. Branch- 
commutativeness is a much weaker condition than commutativeness introduced by Cheb- 
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otarev [2] , since we demand the commutation only of matrices having the same affix k. We 
demonstrate that branch-commutative matrices can be factorized by Moiseyev's technique 
[TTj , so the branch-commutativeness is necessary and sufficient condition for commutative 
factorization. 

In Section 4 we study the Hurd's idea in our terms. Bypass matrices are introduced. 
They are the matrices connecting the values of unknown function Q + on different sheets. 
If all bypass matrices commute, we call matrix G bypass- commutative. The class of 
bypass-commutative matrices is wider than the class of branch-commutative, for example 
all algebraic matrices, to which Hurd's method can be applied are bypass-commutative. 
We show that for any bypass-commutative matrix G a rational factor can be found, 
transforming G into a branch commutative matrix. 



2 Analytical continuation of the factors Q 

2.1 Problem under consideration 

The initial problem of matrix factorization is as follows: 

Problem 1 For a matrix G(k) defined in a narrow strip along the real axis (—e < 
Im [k] < e) find matrices Q + (k), Q~(k) analytical (maybe except some isolated poles), 
continuous, having algebraic growth in the upper (Im [k] > — e) and lower (lm[k] < e) 
half-planes, respectively, and satisfying the equation 

G(k) = Q + (k)Q-(k). (1) 



Algebraic growth hereafter means that there exists a number I, such that all elements 
of corresponding matrices grow at the corresponding half-plane no faster than \k\ l . We 
cannot expect that the elements will grow exactly as some powers of k, since the solutions 
can have logarithmic behaviour. 

We assume that G(k) is an algebraic function, as it happens in many applications. 
Thus, the function itself and the relation ([1]) can be continued from the strip — e < Im [k] < 
e. 

Besides, we assume everywhere that the determinant of G is not equal to zero identi- 
cally. 

2.2 Notations for bypasses 

Let IZq be the Riemann surface of matrix G{k). Below we shall call k an affix of a point 

(k,G(k))en G . 

Let branch points of G(k) have affixes r+ and r~, where Im [r+] > and Im [r~] < 0. 
Make G single-valued on C by performing cuts going from branch points to infinity. The 
cuts can be chosen as 7+ = (r+, +ioo) and 7" = (r~, — ioo). It is important that the cuts 
should not cross the real axis and each other. As a result, the surface TZq becomes split 
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into several sheets. There is a special sheet of IZq, on which equation (JTJ is assumed to 
be valid. Name this sheet a physical sheet. 

Here and below, the structure of Riemann surface is displayed by graphical diagrams. 
Horizontal lines correspond to the sheets, nodes correspond to branch points, and vertical 
lines link sheets, which are connected at a certain branch point. 

In some artificial cases an affix may correspond to several branch points having differ- 
ent orders. Define for each affix its order n^, which is the least common multiple of all 
orders of branch points with corresponding affix. For example, in Fig. [1] a fragment of a 
Riemann surface is shown. Parameter n for the affix fc is equal to 6. 

K 



Figure 1: Order of an affix 

Introduce a notation for the sheets of TZq- Note that later the same notation will 
be used for the sheets of the Riemann surfaces of Q^. Each point of the surface will be 
denoted by (k){w}, where k is an affix, and {w} is a word describing the path, along 
which the argument k should be carried from physical sheet to a selected sheet. The 
structure of this word is explained below. 

Denote bypasses about points rf in positive direction by letters a, and bypasses about 
points r~ in positive direction by bi (Fig. [2]). 




Figure 2: Notation for bypasses 

A series of consecutive bypasses will be denoted by a word of letters a« and 6j. The 
word must be read from left to right, i.e. the first performed bypass corresponds to the 
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left end of the word and the last bypass corresponds to the right end. By default a series 
of bypasses begins from the "physical sheet" . A trivial bypass will be denoted by letter e. 

Define the composition of words w and v as the bypass performed along the way 
composed of w and v. The bypass w is performed first. Denote this composition by wv. 
Let W be the set of all words, and let W a , Wb be the sets of words composed only of the 
letters c^, and only of letters bi, respectively. 

Let G(k){e} be the value of function G(k) on the "physical sheet" . Denote by G(k){w} 
the value of G(k) on the sheet that can be reached by performing the bypass w starting 
from the point (k,G(k){e}). 

The set W can be considered as a group of words, a subject of combinatorial group 
theory. Its generators are the letters Oj, 6j, and the relations have the form 

af = e, bf = e. (2) 

As we shall see below, the same relations are valid for the words describing the Riemann 
surfaces of Q ± . 

Relations ([2]) enable one to determine an inverse element for each w G W without 
introducing new letters for bypasses in negative direction (or without using the symbols 
a^ 1 and b^ 1 ). Using below we assume that for each word w there exists a word w' 1 , 
such that ww' 1 = w~ l w = e. 

Let us demonstrate an example of Riemann surfaces for G and Q + . Take matrix G 
from Daniele's paper [I]: 

/ 1 fci-a(fc) \ 
G{k) = ^) fc2+s(fe) , (3) 

\ fci+s(A) / 

where s(k) = ^k\ — k 2 ] ko, k\ and k% are some complex constants. 

In this case the Riemann surface of G(k) has two sheets and two quadratic branch 
points, namely k = ±fc - Let be Re [ko] > 0. Let letter a denote a bypass about ko, and 
letter b denote a bypass about — ko. 

The scheme for the Riemann surface of G is shown in Fig. |3] a. The upper sheet is 
physical (i.e. it contains the "physical" real axis). 

The scheme of Q + corresponding to this problem is shown in Fig. |3] b. The number 
of sheets is infinite, but all branch points are of second order, and the positive physical 
half-plane contains no branch points. This structure can be revealed, e.g. from jl]. 

2.3 Truncation operators 

Let be w = a%a 2 ■ ■ ■ a n where a, substitutes an arbitrary single letter. Denote by p 
the maximal number, such that the word aia2 . . .a p G W a . Analogically let m be the 
maximal number, such that a±a 2 ■ ■ ■ a m G W&. Obviously, one of this integers is zero, 
since the first letter of the word is either a,j or by 
Define the truncation operators + and ~ by 

w + = a p+ ia p+2 ...oi n , 
w~ = a m+ ia p+2 ...a n . 
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a) 



{e}—f 



{a}={b}-i 



b) n ^ 



{bob}- 



{b}-L 

{ba}- 



7T 



Figure 3: Diagrams of Riemann surfaces for G(k) and Q + (k) 



For example, applying operators + and to the words w = a^b^, v = 6i& 2 «i a 2 we 
obtain 



w 



-- hb 2 , 
aia 2 , 



w = w = aia 2 bib2, 
v + = v = b^a^, 



w 



w 



2.4 Formulae of analytical continuation 

Consider equation ([T|). Both right and left sides of this equation are analytic functions 
in some neighbourhood of the real axis of the physical sheet. Continue Q + and Q~ 
analytically to this domain and, further, onto some Riemann surfaces. Continue also the 
relation ([T]) onto the Riemann surfaces of G, Q + and Q~ . Obviously, the continuation of 
the relation (Til) can be written in the form: 



Q + {k){w}Q-{k){w} = G{k){w). 



(4) 



At this formula (j3J) has sense only for geometrically fixed bypasses. 

Here we are going to find the formulae of analytical continuation for Q ± , i.e. algebraic 
relations connecting Q ± (k){w} with Q ± {k){e}. 

General formulae of analytical continuation can be written in a recursive form as 
follows: 



Theorem 1 Let Q + (k) and Q (k) form a solution of Problem [IJ Then the following 
relations are valid 

Q + {w} = G{w + } G- 1 ^-} Q + {w+-}, (5) 
Q-{ w } = Q-{w-+} G-^w^} G{w-}. (6) 
(A dependence on k is implied for all functions in |3]], 

The proof is rather straightforward and based on the relations following from the 
regularity conditions 

Q + {w} = Q + {w + }, Q-{w} = Q-{w~} (7) 
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literally denoting that Q + is analytical at the points Tj, while Q is analytical at the 
points rj. Let us prove (jSJ) (relation (jSj) is similar). First, according to (j7|), 

Q+M = Q + {«; + }. (8) 

Then, according to (J4]), 

Q+M = Q + {w + } = G{w + }(Q-{w + })-\ (9) 
According to the second relation of (J7|), 

Q+M = Q + {«; + } = G{w + }(Q-{w + })- 1 = G{w + }(Q-{w + -})-\ (10) 
Finally, according to (jlj) 

(Q-^-})- 1 = (G{«; + -})- 1 g + {u; + -} (11) 

and we get (jSJ). 

Note that for any word to their exists some constant c, such that w/ -1 > = e, therefore 
formula © being repeated several times connects g + {u>} with g + {e}. Analogously, 
Q~{w} is connected with g~{e}. The coefficients are always products of known matrices. 

Analytical continuation in the form (JSJ) has been obtained by Hurd [7] for a particular 
case of a single bypass. Hurd's ideas are discussed later in details. 

Using analytical continuation we can investigate the structure of Riemann surface of 
unknown function Q + . For example, the following proposition can be easily proved: 

Proposition 1 Let G(k) be an algebraic matrix, and let the functions Q + (k) and Q~(k) 
form a solution of Problem 1. Then the functions Q + {k) and Q~(k) can be analytically 
continued onto some Riemann surfaces; both functions have branch points only at affixes 
Ti . The order of each branch point is a divisor of corresponding nf . 

A formal proof can be conducted by induction with respect to the length of the word 
w, which is the argument of Q ± (k){w}. 

Generally, solution of Problem CD is not unique: for example the behaviour of different 
solutions at infinity can be different. However, it is easy to prove that all solutions are 
similar up to a meromorphic matrix factor. 

3 Necessary condition for commutative matrix fac- 
torization 

3.1 Necessary condition in the "check-up" form 

Definition 1 Let G(k) be an algebraic matrix, let its branch points have affixes and 
lie aside from the real axis. Let the sets W, W a , and Wb be defined as described above. 
Riemann surface TZq will be called balanced if for any w G W there exist words w a e W a 
and Wb G VV& such that 

G(k){w a } = G(k){w b } = G(k){w}. (12) 
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An example of a balanced Riemann surface is the Riemann surface of scalar function 
+ \[2 + k 2 with an arbitrary choice of the physical sheet. Besides, a surface of an 
arbitrary matrix function, which is a rational combination of k and several square roots 

J rj — k 2 , is balanced. 

An example of a Riemann surface that is not balanced is the surface of the function 
y/i + k + y/-i + k. 

Definition 2 Algebraical matrix G(k) is called branch- commutative, if for any k the val- 
ues of G on different sheets of its Riemann surface commute, i.e. 

[G(k){wi}, G(k){w 2 }} = G(k){wi}G(k){w 2 } - G(k){w 2 }G(k){ Wl } = 0. (13) 

for any different words ui\ and w 2 . 

To illustrate the definition of branch-commutativeness consider a simple example. 
Take a matrix 

/ k 2k s{k) \ 

G(k) = I 2k k -a(fc) , s(k) = Jk 2 - k 2 . (14) 
\ -s{k) s{k) k j 

Let be Im [k ] > 0. There are two letters, a and b, corresponding to bypasses about k 
and — k . The value of G on the sheet {a} is equal to 

/ k 2k -s(k) \ 
G(k){a} = I 2k k s{k) 

\ s(k) -s(k) k J 

Since TZg has two sheets, to check branch-commutativeness one should check only the 
identity 

[G(k){e},G(k){a}} = 0, (15) 

where G(k){e} is defined by (Tl4"|) . Simple calculations show that ffl5|) is fulfilled, so fTl4|) 
is a branch-commutative matrix. 

The necessary condition of commutative factorization is given by the following theo- 
rem: 

Theorem 2 If a matrix G having balanced Riemann surface admits commutative factor- 
ization 

Q + (k)Q-(k) = Q-(k)Q+(k) = G(k), (16) 
then it is a branch- commutative matrix. 

Proof: The formula of analytical continuation (JH]) has been derived for right factoriza- 
tion. One can obtain a similar formula for left factorization G(k) = Q~{k) Q + (k): 

Q+{ w ) = Q + {w + -} G- 1 {w + -}G{w + }. (17) 

Perform the rest of the proof step by step. Here we mark the statements and make 
some comments if the statements are not obvious: 



8 



1. For any word w Q + {w}Q {w} = Q {w}Q + {w} = G{w}. It is an analytical 
continuation of (ITB|) . 

2. For any word w [G{w}, (Q + {w})~ 1 ] = 0. This follows from G{w}(Q + {w})~ 1 = 
Q-{w} = (Q + {w})- 1 G{w}. 

3. For any word w [G{w},Q + {w}] = 0. This can be obtained from the previous 
point by multiplication by G~ l at left and right. 

4. For any v £ W a Q + {e}] = 0. This follows from the previous point and ([7]). 

5. For any word w [G{w},Q + {e}\ = 0. Note that for a matrix with balanced 
Riemann surface for any word w there exists a word v £ W a , such that G{v} = 
G{w}, and Q + {v} = Q + {e}. 

6. For any v\ £ W&, v 2 £ W a [G{viV 2 }, G~ 1 {v 2 }} = 0. This statement can be 
obtained by applying left and right analytical continuation formulae to the word 
Viv 2 and by using the previous point. 

7. For any v x £ W 6 , v 2 £ W a [G{ Vl v 2 }, G{v 2 }} = 0. 

8. The statement of the theorem, by noting that for any w\ and w 2 one can find the 
words Vb £ Wfe and v a £ W a , such that G{vt,v a } = G{wi}, G{v a } = G{w 2 }. 

Theorem |2] is an important result of the paper. Note that since the number of sheets 
of G is finite, the necessary condition can be established by checking a finite number of 
matrix identities. 

3.2 Diagonalization and properties of eigenvectors 

Let an algebraic (not necessarily branch-commutative) matrix G(k) have distinct eigen- 
values almost everywhere (i.e. on the whole complex plane excluding several points). 
Represent this matrix in the form 

G(k) = M(k) diag [Ai, ...,X N ] M(k)- 1 (18) 

Here matrix M(k) consists of vector-columns, which are right eigenvectors of G; Ai . . . Ajv 
are corresponding eigenvalues; iV is dimension of G. Normalize the columns of M by 
making all elements of the first raw of M equal to 1. 

Obviously, for obtaining representation (fl8l) one should first solve the characteristic 
equation for G, and then find a solution of an inhomogeneous linear system for each 
eigenvector. 

Denote Riemann surface of matrix M(k) by 71m- Now we have associated with a 
matrix G two Riemann surfaces: TZq and TZm- Typically, say for Khrapkov matrices, TZm 
has a structure very different from TZq- 

A lot of authors studied matrix factorization problems by formulating a functional 
problem on a Riemann surface. It is important to mention that most of them had in mind 
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the surface TZm, not TZq- Typically, the surface TZm is studied in Moiseyev's context, and 
IZq in Hurd's one. 

Obviously, Riemann surface for the eigenvalues Xj(k) should contain branch points of 
both structures, i.e. of G and of M. 

Let G(k) be a branch- commutative matrix. In this case the set of normalized eigenvec- 
tors must be the same on all sheets of TZm- Therefore, matrix M(k) possesses an important 
property: any bypass about branch points leads to a permutation of the columns, i.e. an 
analytical continuation of each column along a closed contour c on C is some other column 
of M. 

As an example, consider matrix ([3]), which is branch- commutative. As it was men- 
tioned, it has only two branch points, namely ±fc . The scheme of Riemann surface for 
this matrix is shown is Fig. |3] a. It is easy to find that matrix M for this G is as follows: 

( 1 1 \ 

M{k) = ygrg^g vjFggf • (19) 

Matrix M has four branch points, namely ±a/&q — k\ and ± a/ k\ — k\. Generally (i.e. if 
k\ 7^ and ki ^ 0) the branch points of TZm are different from the branch points of TZq. 
The scheme of TZm is shown in Fig. HI 



(k k { ) 



2x1/2 



(k k 2 ) 



( 



kiy 1 



(kl 



2x1/2 



Figure 4: Diagram of TZm 



A transition from one sheet of TZm to another leads to a permutation of the columns 
of M. 



3.3 "Ansatz" form of necessary condition 

Theorem 3 LetG be a branch- commutative matrix NxN, whose eigenvalues are distinct 
almost everywhere. Then it can be represented in the form 

N-l 

G=Y,9m(k)A m (k), (20) 

m=0 

where A(k) is a rational matrix, g m {k) are algebraic functions. Vice versa, any matrix 
admitting a decomposition of the form j[W\) is branch- commutative. 

Proof: The second part of the theorem is obvious, so we are concentrating our efforts on 
the first one. Consider matrix M(k). Let n c be a permutation of columns of M occuring 
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when the argument is carried along a contour c on C starting and terminating at k. Let 
Il c be a matrix containing only numbers and 1, describing permutation n c in matrix 
language, i.e. 

(He)™ = <$ m ,7r c (n), ( 21 ) 

and the permutation of columns of M looks like M — > MU C . 

Construct N functions f m (k), m = 1 ... N as follows. Take N constants /?i . . . (3n such 
that the combinations 

JV 

f m (k) = J2^( M ^ ( 22 ) 

n=l 

almost everywhere obey the relation f mi (k) ^ fm 2 (k) as m i m i- (Here (M)™ are the 
elements of M.) Obviously, f m — >■ /^(m) when the argument is carried along c. 
Construct a combination 

A{k) = M{k) diag [/!(*), . . . , f N (k)} M' l (k). (23) 

Note that the diagonal matrix obeys the relation 

diag [/ We (ij , ■ ■ ■ , U C (N)] = n^ 1 diag [ft,..., f N ] U c . (24) 

Substituting ( |2TT) and ( 124"]) into ( 1231) . conclude that A remains unchanged after any by- 
pass c. Since A is an algebraic matrix by construction, it should be a rational matrix. 

Finally, let us show that G can be expressed in the form (T201) with matrix A constructed 
above. The matrix composed of the elements (F)™ = f™~ x (here m — 1 is a power, 
m = 1 . . . N) has a non-zero determinant almost everywhere. In the opposite case it 
would happen that iV distinct numbers are roots of a polynomial of order smaller than 
N. Therefore any set of iV numbers, for example the eigenvalues of G, can be represented 
as 

N 

X n (k) = Y,9m(k)fr\k) (25) 

m=l 

for almost all k. By construction, g n are algebraic functions. 
The theorem is proved. 

The form (1201) is close to that of [2U|, however on one hand we impose no restrictions 
on the behaviour of the matrices Q ± , and on the other hand, we do not specify the form 
of equation, which matrix A should obey. 

Theorem [3] states that there are two alternative ways to check, whether a diffraction 
matrix G can be factorized commutatively: 1) by checking whether a matrix can be 
represented in a certain form, and 2) by checking conditions (fT3"|) between different sheets. 
The second variant seems more easy. 

Note: Commutative factorization of the matrices having form (|20|) was considered in 
where an explicit formula for the factors was constructed. 
Thus, branch-commutativeness (or, alternatively, the form (|20|) ) is a necessary and 
sufficient condition form commutative factorization of matrices with balanced Riemann 
surfaces. 
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4 Bypass matrices and Hurd's method 



4.1 Bypass matrices 

Let G be an algebraic matrix with a balanced Riemann surface, and let W be the set of 
words associated with this matrix. Let the right factorization problem ([1]) be studied. 

Definition 3 A bypass matrix P w (k) for a word w is defined by the relation: 

P w {k) = Q + (k){w} (Q+(k){e})-\ (26) 

According to the formulae of analytical continuation (jSJ), 

P w (k) = (G{w + }G- 1 {w + -})(G{w + - + }G- 1 {w + - + -}) . . . (27) 

The product in the r.-h.s. is finite, since after some truncations the matrices become equal 
to G{e}. 

Let the number of sheets of TZq be equal to n. Among all bypass matrices we can 
select a finite set of n — 1 basic bypass matrices 

Pj(k) = G(A;){«; j }G'- 1 (A;){e} (28) 

where w± . . . io n _i are any words belonging to Wb, such that all G(k){wj} belong to 
different sheets of TZq, and none of these sheets is the physical one. Any bypass matrix 
can be written as a product of several basic bypass matrices Pj taken in positive or 
negative powers. 

4.2 Hurd's idea and its formalization 

The idea of Hurd [7j can be expressed as follows: sometimes the bypass matrices P w (k) 
can have a structure simplier than that of the matrix G(k). Here we express this simplicity 
in the following form. 

Definition 4 Let G(k) be an algebraic matrix with a balanced Riemann surface. Let P w 
be a set of corresponding bypass matrices. Let all basic bypass matrices commute with 
each other: 

[P j (k),P m (k)] = 0, j,m=l...n-l, (29) 
The matrix G will be called bypass-commutative. 

The approach introduced by Hurd is closely connected with bypass-commutativity. His 
idea was to study the matrix boundary value problem on cuts made from branch points 
to infinity. The matrix coefficient for each cut is one of the bypass matrices. All known 
techniques available for such problems at the current moment require commutation of 
these bypass matrices and their analytical continuations. That is why all known matrices, 
to which Hurd's method was applied successfully are bypass-commutative. 

The class of bypass-commutative matrices is quite wide. For example, all matrices G 
with hyperelliptic Riemann surface are bypass-commutative, since there is a single basic 
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bypass matrix G{a}G 1 {e}. An example of bypass-commutative matrix with a more 
sophisticated Riemann surface is 



where Si = \fk\ — k 2 , S2 = \jk\ — \fk\ — k 2 , k\ and k2 are constants. The Riemann 

surface for such a matrix has four sheets; affixes of branch points are ±k\ and ±a/A^ — k 2 . 
Thus, there are three basic bypass matrices, and their commutativity can be checked 
explicitly. 

Theorem 4 Let matrix G(k) be bypass-commutative. Then there exists a rational matrix 
S(k) such that the matrix G(k)S(k) is branch- commutative. 

Proof. Obviously, if basic bypass matrices commute then all bypass matrices commute. 
Let TZq has n sheets, and the dimension of G be N x N. 
The matrix S{k) can be constructed as follows: 

« 

S(k) = Y,f(k){w,}G-\k){ W] }, (31) 

where Wj is a set of words listing all sheets of TZq- For example, the words from (128]) can 
be taken as Wj for j = 1 . . . n — 1, and w n = e. 

Function f(k) is an arbitrary function, such that it is single-valued on TZq, and the 
r.h-s. of (13TT) has non-zero determinant almost everywhere. As a possible choice, one can 
construct f(k) by the formula 

N N 

/(fc)=/3 ,o + ^^/3 iJ (G'- 1 (fc))} (32) 
i=i j=i 

with almost arbitrary constants 

Since the sum in (13T|) is taken over all sheets of an algebraic function fG, the result 
is single-valued on C, and therefore it is a rational matrix function. 

By construction, the combination G(k)S(k) is a linear combination of some bypass 
matrices or their products. The analytical continuations of G(k)S(k) can also be rep- 
resented as products of bypass matrices. Since the bypass matrices commute with each 
other, G(k)S(k) is a branch- commutative matrix. 

An example (maybe quite simple) of such a consideration can be constructed using 
matrix (I3"UI) . Take function f(k) = 1. Matrix S is constructed by summation of G~ l over 
four sheets of its Riemann surface: 

q(h) = m - k 2 ) fk o\ 

{ } kl + kfk 2 + {k 2 -k 2 )(l- 2k 2 k + 2k*) \ 1 ) 
Simple calculations show that GS is a branch-commutative matrix. 
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Note that Moiseyev's approach can be applied to GS, i.e. an explicit commutative 
factorization can be constructed. 

Note also that if for some matrix G there exists a rational matrix S, such that GS or 
SG is branch-commutative, then obviously all bypass matrices of G commute. That is 
why, bypass-commutativity is a necessary and sufficient condition for a balanced algebraic 
matrix of a possibility to be converted into a branch-commutative matrix. 

5 A short summary 

The main results of this paper are as follows: 

1. Formulae of analytical continuation are derived. 

2. Necessary condition of commutative factorization is found. Namely, a balanced al- 
gebraic matrix should be branch-commutative. This property can be easily checked. 

3. Connection with the "Ansatz" form of the necessary condition is established. 

4. Hurd's method is formalized. It can be applied if a matrix is bypass-commutative. 
It is shown that in this case one can reduce the problem to the commutative factor- 
ization case by multiplication by a rational matrix. 

The necessary condition of commutative matrix factorization for balanced matrix 
is checked as follows. First, one should check, whether the values G(k) taken on 
different sheets commute. If they commute, then the matrix can be factorized by 
Moiseyev's method. Second, one should construct the basic bypass matrices Pj 
and check whether they commute with each other. If they commute, then there 
exists a rational matrix S, multiplication by which transforms matrix G into a 
commutatively factorizable case. 
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